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The evolution of the cosmological perturbation during the oscillatory stage 
of the scalar field is investigated. For the power law potential of the inflaton 
| field, the evolution equation of the Mukhanov's gauge invariant variable is 

o\ '• 

reduced to the Mathieu equation and the density perturbation grows by the 

'. 

parametric resonance. 

O 

Recently, the importance of the parametric resonance is recognized in the reheating 
phase of the inflationary model [1,2]. Due to the coherent oscillation of the background 
inflaton field, the fluctuation of the non-linearly coupled boson field or the fluctuation of self 
interacting inflaton field itself are amplified and the catastrophic particle creation occurs. 
This effect drastically changes the scenario of the reheating considered so far. When one pay 
attention to the evolution of the metric perturbation during the coherent oscillatory phase 
of the scalar field dominated universe, a naive question arises:does the metric perturbation 
undergo a influence of the parametric resonance by the background oscillation of the scalar 
field? 

It is well known that the cosmological perturbation of the scalar field in the oscillatory 
stage has the problematic aspect [3]. If one try to write down the evolution equation for 
the density contrast or the Newtonian potential, the coefficient of the equation becomes 
singular periodically in time due to the background oscillation of the scalar field. This 



behavior is not real. It appears through the reduction of the constrained system to the 
second order differential equation. But it makes difficult to understand the behavior of the 
metric perturbation in the oscillatory phase. 

In this paper, we consider the evolution of the metric perturbation using the gauge 
invariant variable introduced by Mukhanov. The evolution equation for this variable has no 
singular behavior and is suitable to apply the oscillatory phase of the scalar field. We treat 
a spatially flat FRW universe with a minimally coupled scalar field with a potential 
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The background equations are 



2n 

*W = T^(fJ > (n = U,...) (1) 



* 4.4 ( <t> 



H 2 = ^l<P 2 + Vm (2) 
H = (3) 
+ 3H0 + = 0, (4) 

where k = 8nG. The scalar field oscillates if the condition <C m p \ is satisfied. In such a 
situation, we cannot have the exact solution, but using the time averaged potential energy 
and the kinetic energy have the relation 

< 2 >= 2n < V{<j>) >, (5) 

the scale factor and the Hubble parameter can be approximately expressed as 

/ s / t \ 3n n + 1 1 

We define new variables for the scalar field and the time: 

3 / t \ 1 / n 

r] = nt a™+i = nt i^—j , (7) 

<f)(t) = o a"^T0, (8) 

where <f> — 1 at t — to and 0o is a initial value of the scalar field(0o <^ fn P i)- Using the new 
variables, the evolution equation of the scalar field becomes 



w + m 2 n0 2 "- 1 = O, m 2 = ^0 2 . (9) 

For n = 1 (massive scalar field), = cos(m^). For n = 2, = cn{\f\4>Qr]] ^) where cn is a 
elliptic function. 

We use the gauge invariant variables to treat the perturbation whose wavelength is larger 
than the horizon scale. We found that the most convenient variable is 

Q = 50--^ = ^)--^>, (10) 

where tp is the perturbation of the three curvature, 5<j)^ is the gauge invariant variable for 
the scalar field perturbation and \I/ is the gauge invariant Newtonian potential. For the zero 
curvature slice, Q represents the fluctuation of the scalar field. This variable Q was first 
introduced by Mukhanov [4]. As already mentioned, the coefficient of the evolution equation 
for the Newtonian potential or the gauge invariant density contrast becomes singular because 
of the background oscillation of the scalar field. But the evolution equation for Q does not 
have such a singular behavior: 



Q + 3HQ + 



Q = 0. (11) 



The Newtonian potential and the variable Q is connected by the relation 



and the gauge invariant density contrast A which is equal to (^j on the co-moving time 
slice is 

To treat the eq.(ll) more tractable, we change the variable 

3 _ / t \ 1 / n 

Q = a n = nto[ — ) , (14) 



to 



Then 
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Qrm + a "+l 
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Q = 0. (15) 



n + 1 

Using the background equation, we can estimate the time dependence of the each terms in 
this equation: 

( 2-2n> 

.to, 



= n(2n — l)m a "+ 1 



n I I \ I2n-1, 



l-2n> 



-2n N 



As we are considering the situation 77 > to, we can neglect the H 2 , H terms in eq.(15) because 
they are higher order in powers of (1/ 77) compared to other terms. Our basic equation for 
the gauge invariant variable becomes 

1 ~r 



Q-nri 4 



n(2n - l)m 2 2n ~ 2 + k 2 a^ + An(n + l)^ 2n ~Vr, 



Q = 0. 



(16) 



4> is the solution of eq.(9). We first consider n = 1 case(massive scalar). The background 
scalar field is sinusoidal and given by 



Eq.(16) becomes 



m 2 + 



= cos (77777). 



— J sin(m7/) cos(m^) 

a V 



(17) 



Q = 0. 



(18) 



We introduce a dimensionless time variable r = 77177, 

Qtt ~\~ 



1+ f AY_l sin ( 2r ) 
V ma / t 



Q = 0, 

where a oc r 2 / 3 . This equation has the same form as the Mathieu equation: 

Y TT + [A-2qsm(2r)}Y = 0. 



(19) 



(20) 



In our case the coefficient A, q are time dependent functions and the relation between A and 
q is 



(•*)-(!)>. 

Using the stability/instability chart of the Mathieu equation, we can know that the pertur- 
bation will have the effect of parametric resonance of the first unstable band of the Mathieu 
function and grows in time(see figure). We can derive its time evolution by solving eq.(19) 
using multi-time scale method [5]. We introduce a parameter 

6 = i. (22) 

TO 

As the condition tq ^> 1 is equivalent to the condition of coherent oscillation, e is small 
parameter. Rewrite the eq.(19) as 



Qtt ~^ 



1 + 2euJi — e— sin(2T) 

T 



Q = 0, (t>to) (23) 



where lo\ — (—) . We assume the condition (—) < 1 to be the term 2eui small. This 

1 Ze \ ma ) \ ma ) 

means we consider the wavelength larger than the Compton length. We define slow time 
scale r n = e n r. The time derivative with respect to r is replaced by 

jL = D + eD 1 + ---, (24) 

GST 

where D n = We expand 

Q = Q(0) +e Q(l) + .... (25) 

By substituting these expression to eq.(19) and collect the terms of each order of e. From 
the 0(e°) and 0(e 1 ), we have 

O(e ): D 2 o QW + Q^=0, (26) 
0{e l ) : DlQV + q(D = - ^A)AO (0) + 2u*Q<°> - ^ sui(2t)Q< >) . (27) 

The solution of eq. (26) is 

Q(°) = ^(n)e iT + ^*(ri)e- iT . (28) 
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We substitute this to the right hand side of eq. (27) and demand that the secular term which 
is proportional to e lT vanishes. This gives the equation for the amplitude A: 

i^ +Ul A + i^-A* = 0. (29) 
oti At 

Using the definition of e and Ti, we have 

M + U±)\ + U = 0. (30) 
or 2 \ma J t 

Writing A = u + iv(u,v are real), u satisfies the following second order differential equation: 

« = 0. (31) 
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U TT + -Mr + 

3r 



1 / k V 2 



4 V ma I 3r 



Using the eq.(13) and (28), 



_ _ ~~ I ~ ~ u 
A oc Q v (f) v - Q(j) m - -Q v <j> = u + 0(-). (32) 

Therefore the function u is equal to the gauge invariant density contrast A within the 
approximation we are using here. The solution of eq.(31) is 

—- i/e ^((Cfi)- (33> 

where Z u is a Bessel function of order v. We have the critical wavelength Aj = (miJ) -1 / 2 . 
The mode whose wavelength is larger than Aj can grow. If the wavelength is shorter than 
Aj initially, the wavelength is stretched by the cosmic expansion and its exceeds the critical 
length. We can see this behavior by using the chart of Mathieu function. The trajectory 
which started from the stable region moves to the unstable region. For k — > limit, 

Aocr 2 / 3 = a, r" 1 = a" 3 / 2 , 

* oc constant, r" 5/3 = — . (34) 

a 

This behavior is the same as the perturbation in the dust dominated universe. 

For n > 2, the scalar filed oscillation is not sinusoidal. We start searching the solution 
of the equation for y = <f> v : 
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y m + n(2n - l)m z <p 2n ~ 2 y = 0. 



(35) 



Eq.(16) reduces to this equation if the third and the forth term do not exist. We approximate 
the solution of eq.(35) by sinusoidal function:!/ = m sm(cmr]) . c is a some numerical constant 
which defines the period of scalar field oscillation. Using the equation for 0, we have 



n(2n — l)m 



2i2n-2 



y ' 



and using y 



2n(n + l)m 2 2n " V, = (<M>) W , 
msm(cmrj), the equation of Q becomes 



Qrrn "I" 



2 2 4cm . . , 2 4 ("~ 2 ) 

cm sm(2cm?7j + k a "+ 1 



Q = 0. 



(36) 



By introducing the dimensionless time variable r = cmrj, 



Qtt "I" 



4 . s / \ 4 ("~ 2 ) 

1 sm(2r) + — a "+ 1 

r \ cm / 



g = o. 



(37) 



This is also Mathieu equation. It is surprising that this equation contains the n — 1 case if 
we set c = 1. The function A and q are 



4(2-n) 2 

1 + ar 3 ; q = -, 

T 



(38) 



where a = (^) (cnt m) 4( ^ n 2 )/ 3 and we have the relation 

A(q) = ! + «(§)' 



4(2-ra) 
3 



(39) 



Using the chart of Mathieu function, we find that the perturbation also get the effect of the 
parametric resonance of the first unstable band and grows. But as the time goes on, the 
trajectory moves from the unstable region to the stable region and the perturbation will 
oscillate with a constant amplitude. To investigate these behavior, we introduce slow time 
variable and derive the equation for slowly changing amplitude of Q. The procedure is the 
completely same as n = 1 case. The result is 



Q = Ae iT + A*e 

r— + - (—^ 2 

dr 2 \ cm 



4(ra-2) % 

a~^~A+ -A* = 0. 

T 



(40) 
(41) 



The real part of A obeys 

a 2 s,„ o\ /4 2\ 1 
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The solution of this equation is 



. r i(«-2) — I -n — 



3 3J t 



-2 



u = 0. (42) 



-^(sp^r*^)- "- ± 5^V (43) 



For — > limit, we have 



2 /o i\ 2(2n-l) , 3 

Aocra^^aW, r-^a-H+r, (44) 

oc constant, r"3 (4n+1) = — . (45) 

a 

In summary, we found that the evolution equation of the Mukhanov's gauge invariant 
variable in the oscillatory phase of the scalar field can be reduced to the Mathieu equation 
and the evolution of this variable undergoes the effect of the parametric resonance. We can 
interpret the growth of the density perturbation in this phase is caused by the parametric 
resonance. Now we comment on previous works. In paper [6], the analysis is done by using 
the Newtonian approximation which means the wavelength of the perturbation is smaller 
than the horizon length. But the obtained equation for (Sp/p) coincides with the result of 
this paper(eq.(31)). In paper [7], the long wave approximation is used. As the eq.(ll) has 
the exact solution for k — 0, they take in the effect of small k perturbatively and derive 
the evolution of the gauge invariant variables whose wavelength is larger than the horizon 
length. The assumption we used in this paper is the wavelength is larger than the Compton 
length which is smaller than the horizon scale in the oscillatory phase of the scalar field. 
So our treatment is more general. Extension to the non-linearly interacting two scalar field 
system which is a realistic model of the reheating is straightforward and the analysis is now 
going on. We will show the result in a separate publication. 
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FIGURE 

The stability/instability chart of the Mathieu equation (20). The grey is stable and the 
white is unstable region. Three curves shows the time evolution of the parameter A, q for 
the power index n of the scalar field potential V = (jj^ . A = 1 line corresponds to 
k = 0. 
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